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We derive a modified non-perturbative Lorentz-Abraham-Dirac equation, ft satisfies the proper 
conservation laws, particularly, it conserves the generalized momentum, the latter property elim- 
inates the symmetry-breaking runaway solution. The equation allows a consistent calculation of 
the electron current, the radiation effect on the electron momentum, and the radiation itself, for a 
single electron or plasma electrons in strong electromagnetic fields. The equation is applied to a 
simulation of a strong laser pulse interaction with a plasma target. Some analytical solutions are 
also provided. 

PACS numbers: 52.38.-r Laser-plasma interactions, 41.60.-m Radiation by moving charges 
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I. INTRODUCTION 

Lasers now allow us to reach intensities within the focal 
spot of W > 10 22 W/cm 2 Q. Electron motion in fields 
where W 3> 10 18 W/cm 2 , for a typical laser wavelength, 
A ~ 1/im, is ultra-relativistic: 

a = ^4»l, (1) 
mc 

where e is the electron charge, m is its mass, c is the 
speed of light and A is the vector-potential amplitude. 

An accelerated electron in a strong laser field emits 
high-frequency radiation Its back-reaction on the 
electron motion can not be neglected, if in the frame 
where the electron is initially at rest, the energy radi- 
ated during the interaction time is comparable with mc 2 : 
22: fE, 2 c dt > mc 2 , where ctt = J*™ 4 and E is the elec- 
trie field. In the course of a Lorentz transformation, this 
integral transforms proportionally to a wave frequency, 
luq. Indeed, E ~ luqA. The transversal vector potential, 
A, as well as the differential of the wave phase, oj^dt, are 
Lorentz invariant. Since u>q£ — c 2 (ko • p) = c 2 (fco 1 P) is 
invariant, as is the 4-dot-product of the particle momen- 
tum, p l = (f,p), by the wave number, — (^,k ), 
the (sufficient) condition for the radiation reaction sig- 
nificance is as follows: 

ar_ f^cdt = JWdt > mc 2 
mc 2 J 4ir 1.2 kj/(^tm) 2 ~~ £ — cp x 

Here £ and p are the particle energy and momentum 
and the wave is assumed to propagate along the x-axis. 
A high value of the integral in (|2|) may be reached, 
in principle^ at the cost of higher intensity only, W ~ 
10 25 W/cm . In the course of the ELI project (see [H) 
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a laser is expected to reach focusable pulse energy of 1.5 
kJ at A ss 0.8/im, so the radiation effects will be domi- 
nant: JWdt w 2.1 kJ/(^m) 2 > 1.2 kJ/(/im) 2 . Another 
opportunity may be realized while a strong laser pulse 
interacts with energetic electrons, which move oppositely 
to the direction of the pulse propagation. In this case 
£—p x c ss 2£ ^> mc 2 , facilitating the fulfillment of Eq.([2]). 
In such a geometry, powerful X-ray radiation is generated 
in the direction of the electron momentum 3. Ineq.© 
determines the regime, in which the energy is efficiently 
converted to X-ray or 7 bursts. 

In the course of a strong laser pulse interacting with a 
dense plasma the counter-propagating electrons may be 
accelerated in a backward direction by a charge separtion 
field (see, e.g., @). At moderate intensities the genera- 
tion of short pulses of higher-frequency radiation [6| may 
be interpreted as the reflection of the laser pulse from 
these bunches as from a reflecting medium: the frequency 
of the reflected wave, u)( r \ is upshifted: w' r ' ~ (7( m )) 2 u;o, 
if the reflector moves with a Lorentz-factor of ^ 1 
towards an incident wave. Here we consider such high 
laser intensities that emission frequencies are upshifted 
to the hard X-ray and 7 range. At realistic electron den- 
sity, N e < 10 24 cm -3 , the averaged field approximation of 
the reflecting medium is not applicable for emitted pho- 
ton energies exceeding 10 keV, because N e (c/u>^) 3 <C 1. 
The emission from plasma in this case is taken as an in- 
tegral of the radiation intensity from separate electrons, 
rather than as the field of a coherent electric current of 
a 7 range frequency in plasma. Hence, even in a dense 
plasma the emission from separate electrons is essential 
to the analysis. 

While the physical processes involving significant ra- 
diation back-reaction are of growing importance, the 
Lorentz-Abraham-Dirac (LAD) model [?J which is to ac- 
count for this effect, is not free of difficulties, such as 
'runaway' solutions etc. [7|, [8|, [9j, [lfj, [llj, [L2J, [l3|, [l4| . Some 
flaws of the original Dirac version (see Eq.© below), are 
eliminated in the approximate equation, as derived in 
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the book by Landau and Lifshits [g] as Eq.(76.3), see 
also the non-relativistic variant in §75. A slig htly differ- 
ent approximation was found by Eliezer [111 ] , and most 
later versions (see, e.g., [12, E3, EH ) are reducible to those 
listed above. 

New problems arise while considering the transition to 
even higher field intensities, at which Quantum Electro- 
Dynamical (QED) effects come into a power. Delegat- 
ing the discussion of QED strong fields to a forthcom- 
ing publication, we still have to briefly discuss here a 
more general issue of the LAD model conformity with 
the QED principles. Particularly, the nature of gener- 
alized electron momentum (which is substituted for the 
ih\7 operator in the framework of QED) is highly argue- 
able in the LAD model, but this point has hardly been 
discussed so far. Also, there is a controversy between the 
treatment of the radiation in QED as a random process 
of separate photon emission with some probability, and 
the description of the back-reaction of the radiation in 
the LAD model with very smooth functions of time which 
are allowed to be differentiated many times. 

In Section II we describe how we derive the mod- 
ified LAD equation for electrons and account for ra- 
diation from electrons and the electron current in a 
plasma, in a self-consistent manner and in the way which 
does not contradict QED fundamentals. As an applica- 
tion, the basic elements of a particle-in-cell (PIC) nu- 
merical scheme are discussed, which extends the simu- 
lation of laser-plasma interactions to higher intensities 
W > 10 22 W/cm 2 . In Section III we analytically solve 
for the electron motion in a ID wave field in vacuum. 
Results of 3D PIC simulations of laser-plasma interac- 
tion at intensities W ~ 10 22 W/cm 2 are discussed in the 
concluding section. 



II. MODIFIED LAD MODEL 

Here we assume an electron moving in an external elec- 
tromagnetic field and emitting high-frequency radiation. 
In the case of a plasma electron, the external field is the 
averaged field as present in the Maxwell- Vlasov equa- 
tions. We derive the electron momentum equation and 
discuss its (minor) differences from the Dirac theory Q ■ 

Illustration: external field of the ID wave. We 
start from an example of the field of a ID harmonic 
wave with the wave 4- vector, fcg. Recall that the energy- 
momentum exchange of the charged particle with the 
classical field is governed by the Lorentz force, while the 
effect of the emitted/absorbed photons should be inter- 
preted in terms of the photon 4-momentum. The case of 
the ID wave allows both treatments. 

In the course of photon emission with the 4- vector, k\, 
4-momentum is conserved: p\ = p + nhk^ — hk\ , where 
Pq j stand for 4- momenta of the electron before and after 
the emission. In the classical limit of small recoil, the in- 
crement in the electron momentum, 5p l — p\ — p , should 
be small. Therefore, the condition (p-p) = (£/c) 2 — p 2 = 



m 2 c 2 requires that (Sp-p) = and 5p l = fifcq [^^"j —%k\. 
The second term is the 4-momentum transferred to the 
emission and the first term is the gain in 4-momentum 
the electron obtained from the field. The latter for the 
classical external field reduces to the effect of the field 
tensor, F , on the yet unknown current, e5{dx l / dr): 



-F 



dr (k ■ po) 



Here r is the time in the 'Momentarily Comoving' 
Lorentz Frame (MCLF), such that the spatial compo- 
nents of p l Q vanish. In strong fields as in Eq. |T]) , emission 
characteristics are local functions of the wave field (see 
[l5| , §90) . Therefore, the integration in Eq. © reduces to 
a multiplication by At. Expressing F lk — (dA k /dxi) — 
(dA % /dxk), in terms of the 4- vector-potential and using 
the entities, (k ■ A) = 0, (k ■ k ) = uJq/c 2 



can represent 



(fco-Po) Ppi 



k 2 , = one 



and solve Eq.©: 



S( 



dx l 

~d7 



(po-fifci)/i 

m(f L0 ■ /lo)At : 



fho — 



eF?P l o 



(4) 



In the MCLF the current has only spatial components 
and may be expressed in terms of the Lorentz trans- 
formed electric field, Emolf: $ (^) 



p 2 F 2 

e C MCLF 



fiqJlE M CLF 
eB MCLF AT 



"(/lo'/lo)- Hence, the emission is accompa- 
nied by the displacement of the electron along cEmclf- 
Now we average Eqs.© over the emitted photon pa- 
rameters. The averaging (taking a mathematical ex- 
pectation) is done as a weighted integration over du>\ 
with the differential probability of emission per unit of 
time, dW/drduii , the result being multiplied by At, to 
account for the time integration. In the MCLF, aver- 
aging of hcoi gives I At, where I = hwi g^g^r duj\ 
is the total emission intensity. Below we use its ratio 
to the dipole emission intensity, Ie = r o e2 -£'MCLF/ TO ' 

TO 



2e 2 /(3mc 3 ) ~ 6.2 • 10 24 s. In an arbitrary 



frame of reference, averaged Tik\ is the 4-momentum of 
emitted radiation, (^-) ra dAr, expressed in terms of /: 

(^r)rad = (see Q, §73). Analogously, averaging 

(po ■ hki) gives ml At, so that the averaged Eq.(|4]) reads: 

(^r)rad = T oj^^- and the momentum equation for elec- 
tron becomes: 



dp^_ 
~d7 



= eF 



ik Pk 



P l l 



roe 



2 I F* k F klP l 
(mc) 2 



(5) 



where the terms on the right hand side are: the Lorentz 
force, f l L0 , the 4-momentum of the emitted radiation, 

— (^r)rad, and the external field effect, F\,J k ad /c, on the 
current, J t fe ad 



e(%)rad- Multiplying Eq.© by Pi we 



see that d(p ■ p) / dr — 0, maintaining the entity, (p ■ p) 



1 2 

m c . 



General case of an arbitrary external field. 
Eq.© is not specific to the ID wave case and can be 
derived for an arbitrary external electromagnetic field. 
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Seeking the last term in the form of F k J k ad , which is 
mandatory for the 4-momentum exchange with the clas- 
sical external field, and requiring the conservation of (jp-p) 
we obtain Eq.([5]) directly and with no extra assumption. 
Electron current. Now we re- write Eq.© in terms 

of the total electron current, e^fe- = + e ( ) : 

V / 7 



rad 



dp 1 


= -F i 


k dx k 


Ip l 


~d7 


c 


dr 


mc z 


dx l 


= P - + 


I 


eF ik p k 




m 


T °I E 


m 2 c 



(6) 



(7) 



Integrating by volume the equation for the energy- 
momentum tensor for the external field, dT* k t /dx k — 
— \F lk jk, and representing the volume integral, 
J j l dV, of the point-wise current density, j % = 
ecj ^-6 4 (r k — x fc (r))dr, (r k being the coordinate 4- 
vector in an arbitrary Lorentz frame) in terms of e 



we find: £fT^dV = -fi**^ (cf 



dx 1 
dr ' 



§33). 



Eqs.([6][7]) conserve the total energy-momentum. 
V _ 



Hence, 

dp* , 



The generalized momentum may also be conserved. 
However, for a radiating electron this conservation takes 
place, if not only the external field is constant along 
some direction: (n ■ W)A l — 0, but also the projection 

of the emitted momentum, (^r) ra d, onto n vanishes: 
I(p • n) = 0. If the latter condition is not fulfilled and 
I(p • n) ^ 0, then the change in the generalized momen- 
tum, V % — p % + eA % /c, is as follows: 



dV 
dr 



75T =-(»'P) 



(8) 



Eq.© also follows from the quantum relationship, 
n • 5V = —Tin ■ ki, (a conserved generalized momentum 
corresponds to the constant gradient of the electron wave 
function phase along n - see [HI). 

Discussing possible choices of I, we note that the ratio 
I /Ie should be bounded at Ie — > 0. Although to take 
/ = Ie is physically reasonable, there are other interest- 
ing options. Particularly, / can be a random function 
with its average equal to Ie (to trace the quantum the- 
ory limit or to include emission with large photon en- 
ergy). One can apply I expressed in terms of the mod- 
ified emission probability, to treat the processes (like a 
gyrosynchrotron emission, see §90 in |15|) in very strong 
fields, such that the QED effects are not negligible. Af- 
ter all, / can differ from Ie by a choice of w m ; n 7^ 0. 
The latter approach allows us to separate, if desired, the 
high-frequency emission from a lower-frequency averaged 
external field, in simulating laser-plasma interactions. 

Eqs.([5][7|) and their properties result from the asser- 
tion that the electron while emitting moves not strictly 
along the direction of its momentum. Particularly, in the 
MCLF the electron while emitting is not at rest and dis- 
places along eEMCLF- In the MCLF the external electric 



field produces a work at a moving charge, which entirely 
balances the emitted energy: e(^ • E MCLF ) = /. 

The model applicability is limited by the require- 
ment for the current (dx % / dr) r ^ to be essentially non- 
relativistic, which is fulfilled as long as tqI 2 /Ie <C 
mc 2 . To neglect QED effects, the field should be weak: 
eEMCLF <C mc 2 {mc/Ti), and Tols/mc 2 <C (e 2 /Tic) 2 <C 1. 

To compare with the radiation force model we 
use MCLF and approximate within a short time interval 
p = and put / = Ie in spatial components of Eqs. (|6l7p : 
dp/dr 1=3 eE + r e 2 [E x B]/(mc), dx/dr ps (p + r eE)/m. 
Formally, the latter is equivalent to the Newton equation 
with the approximate radiation force as described in 
dill: md 2 x/dr 2 = eE + T (edE/dT + e 2 [E x B]/(mcV). 

Now we compare Eqs.([6][7]) with the LAD equation [3]: 



d 2 x l 
~dr 2 ~ 



eF lk dx k 
mc dr 



d 3 J 



dT 3 



To dx l (d 2 x d 2 x\ 

c 2 rfr" Vdr 2 ' It 2 ) ' ( ' 



(cf. [|, Eqs.(76.1-2)). Re-write Eq.® introducing / 
-mTo(d 2 x k /dT 2 )(d 2 x k /dr 2 ), as in % Eq.(73.4): 



dp^ eF lk dx k 
~dr~ = 



dx l I 



Pb 
m 



dx l 
~dr~ 



t - 



d 2 x i 



dr dr c 2 ' m dr u dr 2 

Comparing this with Eqs.([S][7]) we find that both our 
model and the Dirac theory, as well as the modified ver- 
sion of Eq. ^ as described in [1, [l3|, [l4| (which approx- 
imates ^£ ^£r^T in the right hand side of Eq.©) 
differ from each other with small terms ~ Tq. 

The key distinction, however, is the choice of the 
electron momentum. An interesting survey jl3| shows 
that this choice in the Dirac theory is ambiguous. It 
is problematic too: for p l — mdxi/dr Eq.© conserves 
(p ■ p) — m 2 c 2 , but the generalized momentum is not 
conserved, however symmetric the external field and the 
radiation may be. Particularly, Eq.® allows this elec- 
tron momentum to change in the absence of the external 
field (the runaway solution, see [1, [H, HU), while the 
conservation of the generalized momentum in the MCLF 
would enforce V = p = 0, as long as (n- V)A l = for any 
n and I — 0. With a different choice of the momentum 
(say, pf-> as introduced above) the generalized momentum 
may conserve, but not (p ■ p). So, the distinction of our 
approach from the Dirac model lies in: (1) the incorpo- 
ration of ~ (^r)rad into the relationship between the ve- 
locity, and momentum, p l instead of the "self- force" 

~ ^:(^r)rad into the force equation for dp 1 /dr and (2) 
the use of a different relativistic formulation, providing a 
different set of exact conservation laws. 

Application to the particle-in-cell scheme. 3- 
vector formulation of Eqs.([5][7]) is as simple as: 

^ ={ L + l[SuxB]-^-(6u-f L ), «u + «5u, (10) 



where: / = /r, u = 



Su- 



^ , f L = eE+f[uxB], and 

(11) 



y/m 2 +p 2 /c 
T t L - U(U ■ f L )/c 2 

m 1 + t (u • i L )/(mc 2 ) ' 
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These equations may be applied to plasma electrons in 
order to simulate laser-plasma interactions at high laser 
field intensity. More precisely the equations should 
be solved for 'particles' consisting of a large number 
of electrons, radiating independently and incoherently. 
Their contribution into an averaged electron current is 
e(u + <5u). The spectrum of radiation from relativistic 
electrons is calculated assuming that an angular distri- 
bution is peaked in the direction of the electron momen- 
tum and can be approximated with the <5-function and 

the frequency spectrum, F(r) = K 5 / 3 (r')dr' , 

r = uji/ujc, uj c — |o; r 7 3 , is momentarily close to that 
from circular motion with a rotation frequency, cu r = 
|pxf L |/p 2 : 

-g-At = ffn- 1 g r | - f 

ditduji \ |p| 



7 2 (^u-f L ) 



At. (12) 



The integral of the spectral function is normalized by 
unity, J F(x)dx = 1. The effect of the radiation on the 
'" dfldcdi = 7 2 (<5u 'ft), is entirely 



electron motion, J , n l 
included into Eq.ljlDb. 



III. ELECTRON IN THE ID WAVE 

In the case of the ID wave external field, the electron 
motion can be solved analytically. With the external 
field being a function of £ = (ko ■ x), the relation be- 
tween £ and t is given by a product of Eq.([7]) by ko: 
d£/dr — (ko •p)/m. Multiplying Eq.© by ko, express- 
ing the derivative over r in terms of that over £ and 



assuming I 

d(ko-p)/dt: 
1 



Ie = T~o(k ■ p) 2 c 2 \da/d^\ 2 /m, we obtain: 
-{ko -p) 2 T\da/d£,\ 2 /m, and (cf. to Eq.(2)): 



1 



(ko-p) \{ko-p)Jc = o 



TO 



da(&) 



da 



d£i- (13) 



The transverse momentum, p±, is solved from Eq.(|5]): 



{ko ■ p) 



(ko ■ P) 



-CTo 



da. 



dfi- (14) 



To compare with the Dirac solution for a short pulse 0] , 
consider a single-period symmetric wave: a = a sin(£), 
< £ < 2tt. In the frame of reference in which the 
electron was at rest prior to the interaction the transverse 
components of the electron momentum vanish after the 
pulse: not only they do not turn to infinity, as they would 
in the Dirac runaway solution, but the conservation of 
the generalizied momentum in conjunction with the pulse 

symmetry (J* a 



d£i = 0) entirely eliminates pj_ 
at £ > 2tt. For a pulse of moderate intensity such that 
the left hand side of Eq.(2) is much less than unity, the 
electron gains a small momentum, p x = a? J E 2 dt/(47r), 
in the direction of the pulse propagation, this momentum 
is absorbed from the pulse. The energy, cp x , is absorbed 
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FIG. 1: Energy of backward emitted photons with TuJ > E, 
as a function of E, where a circularly polarized wave of the 
amplitude of ao = 50, interacts with a counter-propagating 
electron of the energy of £ = 180MeV, for different pulse 
durations: 2T (curve 1), 14T (2) and 81T (3), T w 2.7 fs. 



from the pulse and almost equal energy is emitted. For 
strong pulses satisfying Ineq.(2), we present an integral 
spectrum of emission in Fig.l. We see that for a longer 
pulse duration the spectrum is softened. This is a result 
of the radiation reaction: without it the spectrum would 
have the same shape as that of the shorter 5-fs pulse and 
would only increase proportionally to the pulse duration. 



IV. SIMULATIONS AND DISCUSSION 

To demonstrate more realistically the role of the radia- 
tion back-reaction in the laser-plasma interaction we per- 
form a 3D PIC simulation for a 10-cycle linearly polarized 
laser pulse having a step- like profile along the pulse direc- 
tion including 2-A rising and falling edges, and a Gaussian 
profile in the transverse direction, focal diameter 5A, and 
amplitude ao = 70. The laser pulse is incident normally 
on a plasma layer of 10- A length and density no = 3n cr 
where n cr = mc 2 ir/(\e) 2 is the critical density. The sim- 
ulation is performed in the box 20A x 20A x 20A with 
spatial resolution A/20 and 8 electrons per cell, requiring 
in all 6.4 x 10 7 grid cells and 2.6 x 10 s particles. The 
plasma layer is located after a 5- A vacuum layer. Here 
ions are immobile and the time step is At = A/ (40c). 

In this simulation, which corresponds to the intensity 
10 22 W/cm 2 for A = 0.8^m, by the instant t = 20A/c the 
laser pulse loses ~ 27% of its energy, converting ~ 0.9% 
of the incident radiation (or ~ 3.2% of the lost energy) 
to the backward scattered high-frequency radiation. The 
angular distribution of the radiation exceeding 150keV is 
shown in Fig. 2(a). The emitted radiation has a wider an- 
gle along the direction of a transverse electric field. The 
total energy of the pulse equals ~30 J, and emitted back- 
ward high-frequency radiation accounts to ~0.26 J, with 
0.24 J of photon energy above 150 keV (see Fig. 2(b)). 
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FIG. 2: Results of 3D PIC simulation for a linearly polarized 
laser pulse with amplitude ao = 70 entering a soft plasma 
layer (no = 3n cr ). (a) Angular distribution of the backward 
scattered radiation with photon energy above 150 keV and (b) 
total emitted energy as a function of the cut-off frequency. 



The total radiated energy may be close to the particle en- 
ergy for some electrons (see [Hi)- Specifically, this relates 
to a counter-streaming flow of electrons, with momentum 
up to 150mc propagating in the region of the laser pulse, 
as we observe in the simulation. These fast electrons are 
generated in the charge separation field from the flow of 
cold electrons; due to the loss of their energy in the laser 
field, combined with an action of the charge separation 
field of opposite sign, they reverse their motion. However, 
only a minor fraction of electrons counter-propagate with 
high energy, and of these, only a fraction moves in the 
region of the strong laser field, where they can radiate. 
By this account the overall conversion efficiency is dimin- 
ished and does not exceed the order of a few percent. 

We have also applied the described model within a 
particle-in-cell code to simulate processes pertinent to 
fast ignition at laser intensities W > 10 22 W/cm 2 [l7j |. 
The processes in higher fields, such that the corrected 
QED probabilities should be used to simulate the emis- 
sion, are to be considered in a forthcoming publication. 
We also plan to use high contrast pulses from the Her- 
cules laser to drive high-density targets with intensities 
> 10 22 W/cm 2 . We hope that such a study may improve 
understanding of ultra-intense laser-plasma interactions 
and may result in short X- or 7-burst production. 
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